Pseudo-Path Semiclassical Approximation to Transport through Open Quantum 
Billiards: Dyson Equation for Diffractive Scattering 
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We present a semiclassical theory for transport through open billiards of arbitrary convex 
shape that includes diffractively scattered paths at the lead openings. Starting from a Dyson 
equation for the semiclassical Green's function we develop a diagrammatic expansion that allows a 
systematic summation over classical paths and pseudo-paths which consist of classical paths joined 
by diffractive scatterings ("kinks"). This renders the inclusion of an exponentially proliferating 
number of pseudo-path combinations numerically tractable for both regular and chaotic billiards. 
For a circular billiard and the Bunimovich stadium the path sum leads to a good agreement with 
the quantum path length power spectrum up to long path length. Furthermore, we find excellent 
numerical agreement with experimental studies of quantum scattering in microwave billiards where 
pseudo-paths provide a significant contribution. 
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I. INTRODUCTION 

Semiclassical approximations are among the most use- 
ful tools in describing and analyzing ballistic transport 
in mesoscopic systems. On a fundamental level, semi- 
classical techniques allow to build a bridge between clas- 
sical and quantum mechanics: the classical paths carry 
an amplitude which reflects the geometric stability of the 
orbits and a phase that contains the classical action and 
accounts for quantum interference [1-3]. 
Ballistic transport through billiards has been studied ex- 
tensively in the last decade [4-21] and a variety of semi- 
classical approximations [5-13, 22, 23] have been intro- 
duced in order to provide a qualitative and, in part, also 
a quantitative description of these systems. In partic- 
ular, universal conductance fluctuations (UCF) and the 
"weak localization" (WL) have been studied [12] in or- 
der to delineate characteristic differences in the quantum 
transport of classically chaotic and integrable billiards. 
Very recently, quantum shot noise [24] in ballistic cavi- 
ties that are either chaotic [25, 26], regular [27] or dis- 
play a mixed phase space [28] has been used as a probe 
of the quantum-to-classical and chaotic-to regular cross- 
over [25-29]. 

The approach of the (semi) classical limit of ballistic 
quantum transport is both conceptually as well as nu- 
merically non-trivial as it represents, generically, a multi- 
scale problem. The two-dimensional quantum billiard (or 
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quantum dot, see Fig. 1) is characterized by an area A 
or linear dimension D = \f~A~. The quantum wires (or 
leads) to which the billiard is attached have the width 
d. In order to reach sufficiently long dwell times such 
that differences between transiently regular and chaotic 
motion become important, the relation d/D <C 1 should 
hold. To approach the semiclassical limit for the motion 
inside the billiard requires Ajj <C D (Ajj de Broglie wave 
length) or equivalently kD ^> 2tt. Furthermore, if the 
(disorder) potential inside the dot varies over a length 
scale ap, we should require Ac < ap for a semiclassi- 
cal approximation to hold. These conditions pertaining 
to the dot are necessary but not sufficient. Since the 
scattering (S) matrix maps asymptotic scattering states 
onto each other, also the entrance and exit channel states 
in the quantum wire should reach their classical limit, 
Arj *C d or kd 3> 2-7T. The latter limit is virtually impos- 
sible to reach, neither experimentally for quantum dots 
[25, 30, 31] or microwave billiards [32] nor numerically 
[33]. We will therefore focus in the following on the "in- 
terior" or "intermediate" semiclassical regime pertain- 
ing to the interior of the billiard, A^> <C D with con- 
vex hard-walled boundaries such that quantum diffrac- 
tion in the interior can be neglected, with the under- 
standing, however, that quantum effects due to the cou- 
pling to the asymptotic quantum wires have to be taken 
into account. Accordingly, the term "semiclassical ap- 
proximations" refers in the following to approximations 
to the constant-energy Green's function for propagation 
in the interior of the dot G(r, r*,fc) by an approximate 
semiclassical limit, G sc to be discussed below. The 
projections of G onto asymptotic scattering states with 
transverse quantum numbers mT,(mii) yield the ampli- 
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tudes for transmission T mi ,m R from the entrance (left) 
to the exit (right) lead and reflections R mLtTn i . Stan- 
dard semiclassical approximations to G face several fun- 
damental difficulties [8, 9, 12, 16, 23, 34]: among many 
others, unitarity is violated with discrepancies in some 
cases as large as the conductance fluctuations the the- 
ory attempts to describe [9, 10]. Likewise, the anti- 
correlation <5|T| 2 = — <5|-R| 2 between transmission fluctu- 
ations, <5|T| 2 , and the corresponding fluctuations in the 
reflection, <5|i?| 2 , as a function of the wave number k is 
broken. Also, the "weak localization" effect is consider- 
ably overestimated [14, 34]. These difficulties are due to 
the fact that hard- walled billiards possess "sharp edges" 
at the entrance and exit leads even though the interior 
of the dot features a smooth (in the present case, a con- 
stant) potential. At these sharp edges the contacts to the 
quantum wires feature spatial variations of the potential 
where the length scale ap approaches zero. Consequently 
the semiclassical limit Xp> jap <C 1 cannot be reached, no 
matter how small Xd (or large k) is. In other words, the 
quantum properties of the leads influence also the semi- 
classical dynamics in the interior. This observation is 
the starting point for diffractive corrections such as the 
Kirchhoff diffraction [5] or Fraunhofer diffraction [8] . 

We have recently developed a pseudo-path semiclassi- 
cal approximation (PSCA) [13] with pseudo-paths that 
result from spawning of classical paths due to diffractive 
(i.e. non-geometric) reflections in the lead mouths (or 
point contacts). Pseudo-paths play an essential role when 
incorporating indcterministic features into the semiclas- 
sical description of transport. Their existence has also 
recently been pointed out in Ref. [35], although no ex- 
plicit numerical investigations were performed. While 
classical trajectories are either ejected through the exit 
lead contributing to T or return back to the entrance 
lead contributing to R, a quantum wavepacket will do 
both. Pseudo-paths interconnect otherwise disjunct sub- 
sets of classical paths that exit either through the left 
or right lead. The lack of this coupling is responsible 
for violation of the anticorrelation of transmission and 
reflection fluctuations, ST 2 ^ — SR 2 , in standard semi- 
classical approximation (SCA). Likewise, the standard 
SC approximation is expected to fail for quantum shot 
noise [25] that is a signature of this quantum indeter- 
minism. For the special case of the rectangular shaped 
dot with lead openings placed at the midpoints, sum- 
mation of exponentially proliferating pseudo-paths could 
be accomplished by tracing trajectories in an extended 
zone scheme [13, 36]. For arbitrarily shaped billiards 
and, in particular, chaotic billiards where already clas- 
sical trajectories proliferate exponentially, systematic in- 
clusion of paths and pseudo-paths up to the same length 
is considerably more complicated. In the following we 
address this problem within the framework of a semi- 
classical Dyson equation. We present a diagrammatic 
expansion that allows a systematic summation of clas- 
sical path and pseudo-path contributions. Applications 
to the circular billiard (as a prototypical regular system) 




FIG. 1: An open arbitrarily convex-shaped billiard with two 
narrow leads of equal width d: left (L) and right (R). 

and to the Bunimovich stadium billiard (as the prototype 
system for chaotic scattering) show good agreement with 
the numerically calculated exact quantum path length 
spectrum. We furthermore apply the PSCA to recent 
experimental studies of quantum scattering in microwave 
billiards in which pseudo-path contributions could be ex- 
perimentally identified. 

The plan of the paper is as follows: In section II we 
briefly review the standard semiclassical approximation 
and previous attempts to include diffraction effects. The 
pseudo-path semiclassical approximation for arbitrarily 
convex shaped billiards will be presented in section III 
employing a semiclassical version of the Dyson equation. 
We develop a diagrammatic expansion of G in terms of 
paths and pseudo-paths. Its evaluation is numerically 
facilitated by an algebraic matrix representation that re- 
duces path summations up to infinite order to a sequence 
of matrix multiplications and inversions, as discussed in 
section IV. Numerical results and comparison with the 
full quantum results as well as microwave experiments 
are given in section V, followed by a short summary and 
outlook onto future applications in sec. VI. Details on 
the Fraunhofer diffraction approximation are given in the 
appendix. 

II. STANDARD SEMICLASSICAL 
APPROXIMATION 

The conductance of a ballistic two-terminal system (as 
depicted in Fig. 1) is determined by the scattering am- 
plitudes T m through the Landauer formula [37] , 

o 2 / M M \ 
ff(*) = -f E E \ T ^, mL {k)\ 2 , (2.1) 

\m L = l m R = l / 

where M is the number of open modes in the leads (quan- 
tum wires) and T mR ^ mL (k) are the transmission ampli- 
tudes from the m th mode in the entrance lead [referred 
to in the following as left (L) lead] to the (m') th mode 
in the exit lead [referred to as the right (R) lead]. In 
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the following we choose local coordinate systems (xi, yi), 
i = L, R for the leads, where Xi denotes the longitudinal 
and yi the transverse direction of lead i. For simplic- 
ity we use the coding (L,R) for the entrance and exit 
channels throughout this publication, irrespective of the 
actual location at which the leads are attached. The 
projection of the retarded Green's propagator onto the 
transverse wave functions <fi mL (yL) and (l> mR (y R ) of in- 
coming and outgoing modes [12] serves as starting point 
for most scmiclassical theories which approximate trans- 
mission amplitudes from mode to mode m R : 

T mL ,m R 

(k) = -iy/k k,n L ,x L J dy R J dy L <f>* mn (y R ) 

xG(x R ,y R ,x L ,y L ,k)(j) mL (y L ) (2.2) 

Here and in the following we use atomic units (h = |e| = 
m e = 1). The transverse wavefunctions (j) m (y) are given 
for zero magnetic field (B — 0) by 



lKJ> V d \ sm(?fy) to even, 



(2.3) 



where d denotes the width of the leads. We assume, 
for simplicity, that the two leads have identical width. 
The scattering state in the lead is the product of the 
transverse lead wave function 4> m (y) and a plane wave 
in the longitudinal direction with wave vector k XtTn = 
yjk 2 — \mn/d\ 2 . Both leads have a total of (to = 
1,...,M) open (i.e. transmitting) modes. Analogously, 
the amplitude for reflection from the incoming mode to_l 
into a different mode m' L in the same quantum wire is 
given by 



Rm,L,m' L — fimL,m' L * ^Jk XL k mL tXL J d V L I dyL 

4>* m ' L (y'L)G(x' L ,y' L , x L: y Ll k)(i> mL (y L ) . (2.4) 

Equations (2.2) and (2.4) can be written in terms of flux 
normalized projectors onto the left (right) lead or point 
contact (x° L R , y L R ) with matrix elements 

(r\P L \r') = sJ\k XL \ \k' XL \S(f-f')S(x-x° L ) (2.5) 
[e(y-y L + d/2)-Q(y- y ° L -d/2)] 



as 



(r\P R \r') = \KJ 5(r-f')5(x-x R ) (2.6) 

[e(y-y R + d/2)-e(y-y R -d/2)] 

T mR ,m L {k) = -i{m R \P R G{k)P L \m L ) (2.7a) 
i(m' L \P L G(k)P L \m L ). (2.7b) 



The term \k XL \ {\k' x \) denotes the longitudinal compo- 
nent of the momentum of the incoming (outgoing) wave 
function at the left lead. 

The semiclassical approximation to the scattering am- 
plitudes is obtained by approximating the Green's func- 
tion G (r-2, r*i, k) in Eq. (2.2) by the semiclassical Green's 
propagator G sc (r 2 , f\ , k) . The standard scmiclassical 
Green's propagator G sc (r 2 7 ri 7 k), the Fourier-Laplace 
transform of the van Vleck propagator evaluated in sta- 
tionary phase approximation (SPA), describes the prob- 
ability amplitude for propagation from f\ to r 2 at a fixed 
energy, E = k 2 /2. It can be expressed in terms of a 
sum over all classical paths of energy E (or, equivalcntly, 
wavevector k) connecting these two points [1] , 



G sc (r 2 ,n,k) 



sc 



\D a (f 2 ,fi, k)\ 



1/2 



(2.8) 



? - (2tu)V2 

a: n— *r2 

x exp ^ {s a {r 2 ,f u k) - ) 

Here, S a {r2,T\,k) = kL a is the action of the path a of 
length L a . D a (f 2 ,fi,k) is the classical deflection fac- 
tor [1] which describes the stability of the paths and /i a 
denotes the Maslov index of the path a. In line with 
the semiclassical approximation, the double integrals in 
Eqs. (2.2) and (2.4) are frequently evaluated in station- 
ary phase approximation. Physically this means that the 
paths are entering and exiting the cavity only with the 
discrete angles 9 m = arcsin[TO7r/((ifc)] due to the quan- 
tization of the transverse momentum in the leads. For 
completeness we mention at this point that the path-sum 
in Eq. (2.8) contains also those orbit pairs described in 
Ref. [23] which yield a weak-localization correction be- 
yond the diagonal approximation. 

Several strategies have been proposed to introduce 
diffractive effects in order to quantitatively improve the 
semiclassical theory for transport through open quantum 
billiards. A straightforward way is to eliminate the SPA 
for the double integral in Eq. (2.4). Expanding the ac- 
tion in the semiclassical Green's function Eq. (2.8) to 
first order in the transverse coordinate, the integral takes 
the form of a Fraunhofer diffraction integral and can 
be evaluated analytically [8]. On this level of approxi- 
mation, diffraction effects are thereby automatically in- 
cluded, however only upon entering and exiting the cav- 
ity, not during propagation inside the cavity. Schwieters 
et al. [5] employed Kirchhoff diffraction theory to calcu- 
late the diffractional weight of paths entering and exit- 
ing the billiard. In addition, they introduced the concept 
of "ghost-paths": paths that are specularly reflected at 
the lead opening due to diffractive effects. The proper 
use of a diffractions weight for classical paths allowed 
the quantitative determination of the peak-heights in the 
power spectra of the transmission and reflection ampli- 
tudes [5, 8] - at least for short path lengths. Ghost orbits 
could account for some of the peaks that were missing 
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in the semiclassical spectra [5]. Several deficiencies re- 
mained, however, unresolved: unitarity of the semiclassi- 
cal S matrix is, typically, violated; the weak localization 
peak is significantly underestimated, and the semiclassi- 
cal pathlcngth (£) spectrum [15], 



jSC 



rlk p ikl T sc 
un, o ■ L m',m 



(k) 



(2.9) 



fails to account for all the peaks and overestimates the 
corresponding quantum pathlength spectrum for large t. 
For the special case of a rectangular billiard we could re- 
cently demonstrate [13] that a systematic inclusion of 
pseudo-paths within the pseudo-path semiclassical ap- 
proximation (PSCA) has the potential to overcome these 
deficiencies. In the following we derive the PSCA for 
an arbitrary convex billiard from a semiclassical Dyson 
equation and investigate its properties numerically. 




(b) 




— /~> sc p /-» sc p s* sc 

— Ulr ~r Urr ~r Url 



= 



III. DYSON EQUATION FOR THE PSCA 

Starting point is a semiclassical version of the Dyson 
equation for the Green's function where the standard 
semiclassical Green's function G sc plays the role of the 
unperturbed Green's function and the diffractive scat- 
terings at the lead openings (or point contacts) are the 
perturbation. Accordingly, we have 



G 



PSC 



= G 



sc 



g sc vg psc 



(3.1) 



where the perturbation "potential" V is given in terms 
of the projectors Eqs. (2.5, 2.6) as 



V = Pr 



Pn 



(3.2) 



Iterative solution by summation 



G 



PSC 



G 



sc 



i=0 



( VG SCy = G 



SC E 

i=0 



((Pl + Pr)G 



SC\ 



(3.3) 

includes diffractive scatterings into the G PSC to all or- 
ders. 

The key to the multiple diffractive scattering expansion 
is that within the semiclassical expansion each projection 
operator onto the L and R point contacts selects classical 
trajectories emanating from or ending up at the leads 
and at the same time spawns new generations of classical 
trajectories. 

Noting that G PSC will only be evaluated in the domain 
of Pl or Pr, Eq. (3.3) can be reorganized in terms of a 
2x2 matrix Dyson equation. We decompose G sc as 
follows 



Pl G sc P l 



r<sc 

^LL 



E G « 



sc 



Pl G sc P r => 



r*sc 

U LR 



E^ 



sc 



(3.4a) 



(3.4b) 



FIG. 2: Diagramatic expansion of paths. Double line: GfF, 
single line: contribution of individual classical trajectories 
G^f., dots: replica of point contacts along the propagation 

direction k ("time"), (a) Typical term appearing in the ex- 
pansion of G PSC . (b) Expansion of Gff in terms of classical 
path contributions. 



P R G sc P l 



n SC 

^RL 



E G «HI 



(3.4c) 



P R g sc p r 



G S r°r = E GS o° RR (3-4d) 



In Eq. (3.4) the index refers to paths that emanate 
from the point contact j (j = L, R) and end up at point 
contact i(i = L,R). Each of the four disjunct subsets of 
classical paths is, in general, infinite. In the following we 
denote the truncated number of trajectories of the cor- 
responding class by iVy. Each Gf^ containing a large 
finite number Njj or an infinite number of trajectories is 
diagramatically represented by a double line, each contri- 
bution of an individual trajectory by a single line (Fig. 2). 
The trajectories emanate or end on vertices representing 
the L or R point contacts. 

In this formalism Eq. (3.1) now becomes 



r<psc r*PSC 

^LL ^LR 

r<psc r<psc 

^RL ^RR 



I ^LL U LR \ , / ^LL U LR \ 

\ r<sc rsc ) + \ rsc r<sc j 

\^RL ^RR/ K^RL ^RR/ 

(Pl 0\ (G P d C G p LR \°\ 
1.0 Pr) \G%° G™>) 



(3.5) 



Its solution G PSC denoted by solid lines can now be 
diagramatically represented (Fig. 3) as a sum over all 
pseudo-paths that result from the couplings of classical- 
path Green's functions by successive diffractive scatter- 
ings at point contacts. Only pseudopath combinations 
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(or diagrams) contribute that are connected at vertices 
L or R. 

Equation (3.5) can be formally solved by matrix inver- 
sion, 



(G p d c Gl%>\ _ 

\ r<psc r<psc ) — 
V^rl u rr J 



1 - Gf c L P L 
-G S r C lPl 

I <^LL ^LR \ 

\ r<sc rsc j 

\ U RL U RR/ 



- G LR P R 
1 - G RR P R 



(3.6) 



resulting in a sum over (up to) infinitely long pseudo 
paths with (up to) an infinite number of diffractive scat- 
terings. Note that the number of classical paths between 
two diffractive "kinks" and their lengths may reach in- 
finity as well. The handling of this double limit plays an 
important role in the numerical implementation as de- 
structive interference of pseudo-paths and classical paths 
of comparable length must be properly taken into ac- 
count. The evaluation of Eq. (3.5) and (3.6) is not 
straight-forward as each operator product in Eq. (3.3) 
contains a multi-dimensional integral over M 2 . For il- 



lustrative purposes we explicitly give the first-correction 
term (G SC P L G SC + G SC P R G SC ) of Eq. (3.5). We have 
to calculate the double integral e.g., 



r<sc r> n sc 



u; r n ,, l;L = j d 2 r' I d 2 rGi%(? 2 ,?')P R (r',n 

R 2 M 

SC 



x G R Ur,n), 



(3.7) 



where r 2 (fi) are both located on the left point contact. 
The double integral reduces to a one-dimensional integral 
along the lead openings due to the ^-functions in the 
projector Eq. (2.6), 



r<sc n ^sc 

U LR^R U RL 



,1/2 



d/2 



d V \J\ k x R \ \k' XR \ 



(3.8) 



x Gi%(r 2 ,x R , y R + y) G&(a&, Vr + V, 

Inserting the expression for the semiclassical Green's 
function Eq. (2.8), and using the abbreviation r R — 
(x R , y R + y) we obtain 



n sc T3 r<sc 
u lr^R- u rl 



d/2 

EE / *V V& \K R \G s a c LR {f 2 ,r R ) Gf^(rfl,rl) = -^EEP«-I \ D ^\ 



-d/2 



d/2 



dy J\k XR \ \k' XR \exp ik(L ctLR (f 1 ,f R ) + L aRL (r R ,r 2 )) - i- (fi aRL + fi aLR ) 



-d/2 



(3.9) 



where we have assumed that the classical deflection fac- 
tors D a are smooth functions over the range of the lead 
mouth and can be approximated by their value D a at 
the center. In Eq. (3.9) a R L denotes paths connecting 
the left with the right lead (L — > R) while oll R repre- 
sents paths (R — > L). In order to solve the integral in 
Eq. (3.9) analytically, we expand the path length L a of 
the classical paths a leading from point f*j to the lead 
mouth f R (here we set y R — for simplicity and only 



L R 




FIG. 3: Thick solid line: G PSC ; first terms of the expansion 
of the Dyson equation for Gjjf c . 



consider regular billiards): 
L aRL {r R ,n) = L aRL + sm(9 aRL )y-\ y . . . , 

Ll a RL 

(3-10) 

where L aRL = L aRL (f R , fi) is the length of the path that 
reaches the center of the lead mouth. 

Taking into account only the first-order correction in 
the lateral displacement, the integral in Eq. (3.9) takes 
the form of a Fraunhofer diffraction integral and can be 
solved analytically [8, 13] (see also appendix), 

d/2 

I = J dyy/\k XR (y)\ \k' XR {y)\ 

-d/2 

x exp{ik[L aLR (y) + L aRL (y)}} (3.11) 

d/2 

~ fcy / cos(6» Ql/ii )cos(6»^ RL ) J dyexp{ik[L aLR + L aRL ] 

-d/2 

+ ik [sin(0 a£R )+sm(0' R j] y) 
= exp[ik(L aLR + L aRL )]r(8 aLR ,e' aRL ,k), 
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where 9' a is the ending angle of the incoming path and 
6 aLR is the starting angle of the exiting path. Further- 
more, r(6 aLR ,6' aRL ,k) corresponds to the Fraunhofer re- 
flection coefficient at an open lead [13]: 

r{O aLR ,0' aRL ,k) = 2^/cos(0 a£ Jcos 

( r kd /„• a 



X 



sin [!f(smO aLR -+ 

cin -U cr 



sm6 aLR +sin6^ Hi 




Finally, this first order correction term in Eq. (3.7) can 
be written explicitly as 



G LRPR G RL = 2^2^ G a LR r { aLRi O ' aRL i k ) G , 



t cc R l 



(3.13) 

Note that classical paths contributing to G sc do not 
incorporate the finite width of the point contact since 
we take the limit d — > in the simulation of classi- 
cal paths. Instead they are specularly reflected at the 
open lead mouth R as if there were a hard wall. Such 
a path will interfere with the pseudo-path that has al- 
most the same length and topology but experiences a 
non-specular, i.e. diffractive reflection at the same point 
contact. 

All higher-order corrections are evaluated analogously. 
Each additional vertex connecting two Gsc gives rise to 
an additional Fraunhofer integral with an interior reflec- 
tion amplitude r(9',9). 



IV. MATRIX REPRESENTATION OF 
PSEUDO-PATH SUM 

A numerical representation of the Dyson equation 
[Eqs. (3.5,3.6)] requires the truncation of the number of 
contributing paths in each of the four basic Green's func- 
tions [Eq. (3.4)] to large but finite numbers. Let Nll 
denote the number of paths leading from f® L back to r^. 
Accordingly, N RR is the number of paths of the class 
(R — > R), Njil the number for (L — > R) and Nl R the 
number for (R — > L). Note that N R l — Nl R for systems 
with time reversal symmetry. The total number of con- 
tributing paths is N = N LL + N RR + N LR + N LR . In 
order to numerically solve the semiclassical Dyson Equa- 
tion (3.5), we write the semiclassical Green's function 
G sc 

MS el diagonal N x N matrix where each diagonal 
matrix element represents the contribution of one par- 
ticular classical path. Distinguishing the four different 
subclasses of paths, the matrix can be written in the fol- 
lowing form: 



o \ 



G sc (k) = 



v ° 



(4.1) 



where the matrix elements of the sub-matrices are given 
by, e.g., 



G sc (k) 



\D aL 



1 1/2 



{2m) 



1/2 



■ exp 



'( 



kL n 



(4.2) 

and a LL = 1, . . .,N LL . 

Due to diffractive coupling at the lead mouths, the ma- 
trix G PSC that represents the Green's function in the 
pseudo-path semiclassical approximation contains also 
non-diagonal matrix elements. A specific matrix element 
(a 1 , a) consists of a sum of all path-combinations that 
contain the classical path a as first segment and the clas- 
sical paths a' as last segment. 

The coupling between different classical paths is repre- 
sented by the non-diagonal vertex matrix r (k) . The ma- 
trix elements are given - within the Fraunhofer diffraction 
approximation - by the reflection amplitudes [Eq. (3.12)]. 
Since r{6, 9' , k) depends on both the angle of emission, 6, 
of the new path and the angle of incidence, 6', of the 
previous path, 16 block matrices would result. How- 
ever due to the restriction imposed by the projections 
Pl and P R that the endpoint (L or R) of the incom- 
ing trajectory must agree with the starting point of the 
outgoing trajectory, effectively only 8 block matrices can 
contribute. These can be characterized by the starting 
points and endpoints of classical trajectories at the ver- 
tex. The block matrices resulting from the projection Pl 
are ^ LL LL , L LL LR , L RL LL and g RL LR - Analogously, 
from the vertex P Rl we get r 



and lLR,RL- =LL,LR 

the elements 



T T 
= RR, RL —LR. RR 



■RR, RR ' 

is a N LL x N LR matrix with 



L arr . a , (k) = {r(e aLL ,0' aLR ,k)}. (4.3) 



The full matrix is 



r(k) 



r =RL,LL^ 



(k)\ 



=LL,LR 
^RL,LR^ 



V 



^rr,rlW L RRRR {k) 
L LR , RL ( k ) L m (k) 



/ 

(4.4) 

For billiards of arbitrary shape and positions of the leads, 
no further reductions are possible. Only for structures 
with discrete geometric symmetries [L <-> R) or time- 
reversal symmetry, the number of non-equivalent trajec- 
tories and thus of independent amplitudes r{9 1 9 1 ,k) is 
reduced. 

With Eqs. (4.3) and (4.4), the Dyson equation Eq. (3.5) 
can now be written as an algebraic matrix equation 

oo 

Qf sc (k) = Q: sc (k) [L(k)Q: SC (k)] 1 (4.5a) 



1=0 



■g SC {k) [l-L{k)Q SC '(k)\ 



(4.5b) 
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Equation (4.5b) represents the "exact" summation over 
pseudo-paths with up to an infinite number of diffrac- 
tive scatterings. The accuracy of the result is, however, 
limited by the fact that we can only take into account a 
finite number TV of classical paths. 

Finally, calculation of the S matrix elements re- 
quires the projection of G P c onto the asymptotic 
scattering states in the left and right quantum wire 
(Eqs. (2.7,2.8). Following the same line of reasoning as in 
Eqs. (3.9,3.10,3.11,3.12), the projections Pl and Pr give 
rise to a transmission amplitude in Fraunhofer diffraction 
approximation. The transmission amplitude from incom- 
ing mode to to a classical path a inside the cavity with 
launching angle 9 a is given by [13]. 



+ 



2 cos 9 a 
kd 



sin [(fcsin(9 Q + ?f) 



sin 



sin# Q — 



m-ir 
kd 



(4.6) 



Likewise, the transmission amplitude for a trajectory 
a approaching the point contact with angle 9 a to exit in 
mode to is also given by Eq. (4.6). 

With Eq. (4.6) we form now amplitude matrices to map 
the asymptotic scattering states onto the N x N repre- 
sentation of G PSC . The N x 2M matrix for the incoming 
state, where M is the number of open transverse modes 
(m L . R = 1, . . . , M), is given by 



A ln {k) 



tm L {9a LL ,k) 
trriL {@ctRL j k) 

tm R {9a Lr n k) 



(4.7) 



The corresponding projection amplitude for the outgoing 
scattering state is a 2M x 2N matrix and reads 



A out {k) 



tmL^ ctLL ik) imL^fiaLRik) 

t mR (9 aRL , k) t mR (9 aRR ,k) 



(4- 



r 



The 2M x 2M dimensional S matrix follows now from 



S PSC {k) = A out {k)G 1JbU (k)A m {k) 



PSCi 



(4.9) 



Pm' L ,mL ^-ra'^.rafi 
Pm R ,mL Rm' R ,mR 



With Eq. (4.9) the semiclassical calculation of the S ma- 
trix of a hard-walled quantum billiard with two point 
contacts is reduced to a sequence of matrix multiplica- 
tions. System-specific input are the data (length, angle 
of incidence and emission from point contact and Maslov 
index) of each of the four classes (L — > L, L — > R, R — > 
L, R — * L) of classical trajectories. 



V. NUMERICAL RESULTS 

In the following we will apply the pseudo-path semi- 
classical approximation to the calculation of transmission 
and reflection amplitudes in different regular and chaotic 
structures (i.e. circle, rectangle and stadium). We com- 
pare the PSCA calculations with the results of the stan- 
dard semiclassical approximation [8] , exact quantum cal- 
culations [38] and microwave experiments [39] . 

In order to evaluate the scattering matrix elements 



T 



and Rn 



numerically we have used the trun- 



maximum number of diffractive scatterings (or "kinks" ) . 
We therefore relate K to the maximum path-length £ max 
of classical paths included. In order to include all pseudo- 
paths with at least the same length, we require that the 
length of the shortest trajectory in the system l m i n mul- 
tiplied by K + 1 has to be larger than £ max to guarantee 
that all pseudo-paths up to £ ma x are included. Hence we 
require 



K > 



1, 



(5.1) 



cated form of Eq. (4.5a) with i < K, where K is the 



where the bracket stands for the largest integer less 
than the argument. This requirement only assure that 
each classical path is shadowed by a pseudo-path of 
comparable length. The converse is evidently not the 
case. Pseudo-paths are permitted with path-length up 
to £ m ax(K + 1) for which no classical counterpart of 
comparable length is included. We account for this 
deficiency by Fourier filtering the power spectra P(£) of 
the S- matrix elements for £ > £ max . The classical input 
data are determined analytically for the square and circle 
billiards while for the chaotic stadium billiards the data 
are generated by calculating the classical trajectories 
numerically. 

A Regular structures 

For regular structures classical and pseudo-paths can be 
easily enumerated and calculations up to very high path 
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FIG. 4: Power spectra IT22 (^) | 2 of the transmission am- 
plitude in the circle billiard with perpendicular leads with 
R = sj\ + 4/-7T and d = 0.25 (see inset Fig. 5) for a finite 
window of k, 1 < k < 6 in units of n/d. 
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FIG. 5: Power spectra J-/?22(^)| 2 of the transmission am- 
plitude in the circle billiard with perpendicular leads with 
R — yl+4/7T and d — 0.25 for a finite window of k, 
1 < k < 6 in units of n/d. For more details see text. 
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FIG. 6: Power spectra of the transmission amplitude (m = 
n — 1) in the stadium billiard with perpendicular leads with 
R — 1, I = 2R, and d — 0.25 for a finite window of k, < fc < 
5 in units of n/d. 



lack of individual peaks (i.e. pseudo-paths) present in the 
quantum spectrum but missing in the SCA spectra (e.g. 
at I « 27.5; see white arrows in Fig. 4) can be addressed 
in the transmission amplitude. 

This defect of the SCA can be seen more clearly in 
the plot for R22 (see arrows in the upper part of Fig. 5). 
The pseudo-path semiclassics provides for the additional 
peaks in agreement with the quantum calculation. The 
insets in the lower part of Fig. 5 show schematically the 
geometry of dominant pseudo-paths, which experience 
at least one diffractive reflection at an open lead. The 
number next to the insets gives the information on how 
often the segment of a classical path has to be traversed 
to build up the corresponding pseudo-path. 



length (£ = 40) have been performed as a sensitive test of 
semiclassical approximations. We present in the follow- 
ing results for the path-length power spectrum [Eq. (2.9)] . 
Figs. 4 and 5 show the power spectra of the transmis- 
sion and reflection amplitudes for the second mode, T22 
and R22, in the circular billiard with perpendicular leads. 
Note that all the scattering geometries which we investi- 
gate in the following have the same cavity area A = 4 + 7r, 
or, accordingly, a linear dimension D — y/A + it « 2.7. 

Here we mainly concentrate on the improvement of the 
PSCA in comparison to the SCA. For T22, the overesti- 
mation of long paths (£ > 18) in the SCA calculation is 
clearly visible: compare the peaks labeled by black ar- 
rows in Fig. 4a with those in Fig. 4b. The contributions 
of the additional pseudo-paths included in the PSCA are 
responsible for the cancellation effects which suppress the 
peak heights. The peak heights of the PSCA agree well 
with the quantum calculations for long paths. Also the 



B Chaotic structures (the Bunimovich stadium) 

Up to now we have compared the SCA and PSCA 
for open billiards with regular classical dynamics. We 
turn now to chaotic (i.e. non-regular) billiards. The 
Bunimovich stadium [33, 40, 41] serves as prototype 
system for structures with chaotic classical dynamics. 
Before we discuss our numerical results, we point 
out some of the characteristic differences between the 
regular and chaotic structures. In the stadium billiard 
the number of path bundles [8] up to a fixed length 
increases exponentially. The exponential proliferation 
represents a major challenge for the summation of paths. 
This leads to technical problems for the calculation of 
transport quantities in chaotic structures [33] due to the 
limitation of the number of paths that realistically can 
be taken into account. Furthermore, the path length 
distribution in chaotic structures differs qualitatively 
from that of regular structures. For the latter we find an 
algebraic decay for the classical path-length distribution 
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in contrast to the exponential decay for classically 
chaotic structures. This is to be distinguished from the 
exponential decay behavior for pseudo-paths. Wirtz et 
al. [13] have shown that for the rectangular billiard 
the SCA leads to a linear (algebraic) decay of the path 
length power spectra, while the PSCA gives rise to 
an exponential decay when pseudo-paths are included. 
However, in chaotic systems, already classical paths 
proliferate exponentially as a function of the path length 
and can account for exponential suppression of large path 
lengths. Therefore, the lack of pseudo-paths which also 
proliferate exponentially is less dramatically felt than in 
regular systems where the exponential proliferation of 
pseudo-paths competes with only linear proliferation of 
classical paths. This observation is key to the surprising 
findings in previous semiclassical calculations that the 
agreement between the semiclassical and the quantum 
pathlength spectrum is better for chaotic rather than 
for regular systems [8, 12, 15, 20, 33]. However, also in 
chaotic systems, the effects of diffractive pseudo-paths 
can be clearly seen. Fig. 6 shows the power spectra of 
transmission amplitudes for the stadium with perpendic- 
ular leads. For Tn, remarkably, the path-length spectra 
within PSCA displays fewer pronounced peaks than the 
SCA in agreement with the quantum calculation. The 
reason is that in a chaotic system the high density of 
pseudo-paths effectively causes path shadowing of true 
classical paths. As a result, some of the classical peaks 
are drastically reduced by destructive interference even 
for comparatively short path length. 

C Comparison with microwave experiments 

As a third application we discuss the comparison between 
experimental studies of geometry-specific quantum scat- 
tering in microwave billiards [39] and semiclassical ap- 
proximations (SCA and PSCA). The physics and mod- 
eling of microwave cavities are conceptually similar to 
that of semiconductor quantum dots due to the equiva- 
lence of the time-independent Schrodingcr and Helmholtz 
equations [32]. Moreover, for microwave frequencies 
v < v max = c/2h [32], where h is the height of the mi- 
crowave billiard, only a single transverse mode is sup- 
ported by the cavity. This reduces the electromagnetic 
boundary conditions to Dirichlet-boundary conditions al- 
lowing for an exact correspondence between electrody- 
namics and quantum mechanics. Accordingly, the com- 
ponent of the electrical field perpendicular to the plane 
of the microwave billiard corresponds to the quantum 
mechanical wave function. 

A direct measurement of transmission and reflection 
amplitudes (i.e. electrical field amplitudes) becomes pos- 
sible. Moreover the semiclassical analysis presented 
above is directly applicable to describe scattering in mi- 
crowave billiards. The geometry of the microwave res- 
onator used in the experiment [39] is shown as inset in 
Fig. 7. The transmission and reflection amplitude were 
taken in the frequency range 13 GHz < v < 18 GHz, 
where only one mode is propagating in the waveguides. 




FIG. 7: Power spectra |Tn(^)| of the experimental and cal- 
culated (QM, SCA and PSCA) transmission amplitude for 
the rectangular billiard (L — 225mm, D = 237mm and 
d = 15.8mm) with opposite leads (not centered) for a finite 
window of u, 13 GHz < v < 18 GHz. The characteristic peaks 
are identified in terms of classical transmitted trajectories (in- 
sets). Experimental data by H. Schanze [39]. 
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FIG. 8: Same as Fig. 7, but for reflection. The insets in 
Fig. (b) show the geometry of the pseudo-paths. Numbers 
next to insets give the number of revolutions inside the bil- 
liard. 



Figs. 7 and 8 show the experimental and calculated data 
for the power spectra of the transmission and reflection 
amplitude, |Tii(£)| 2 and |^n(^)| 2 . The important role 
of pseudo-paths in reproducing the correct peak pattern 
in good agreement with experimental results can be seen 
in both the transmission and the reflection spectra. Es- 
pecially in |i?n(^)| 2 distinct pseudo-paths (see insets in 
Fig 8b) appear which were noted in Ref. [39]. For ex- 
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ample, the peak at I « 9 (white arrow in Fig. 8a) is 
caused by the trajectory with the length i sa 4.5 which, 
after one revolution in the billiard, is reflected back at 
the exit by the lead mouth, so that it continues for one 
more revolution. Thus its total length t « 2 x 4.5 = 9. 
Of course, such non-classical trajectories are not included 
in the standard semiclassical approximation (Fig 8a) . In 
this special case, the reflection at the open lead mouth is 
specular or geometric. This peak would therefore also be 
present in the semiclassical approximation suggested by 
Schwieters et al. [5]. By contrast, non-geometric reflec- 
tions ("kinks") by diffractive scattering are also present, 
e.g. the peak at I — 7.5 (see inset in Fig. 8b). The latter 
class is only contained in the PSCA. 



VI. CONCLUSIONS 




5 10 15 20 25 30 ' 5 10 15 20 25 30 

k in units of n/d 



We have presented an extension of the pseudo-path 
semiclassical approximation (PSCA) [13] to billiards 
with arbitrary convex shape. A diagrammatic expansion 
of the semiclassical Dyson equation for the Green's 
function in terms of diffractive scattering diagrams is 
developed. The unperturbed Green's function represents 
the sum over classical paths. Each encounter with the 
lead (point contact) spawns new trajectories. By joining 
disjunct classical paths due to non-geometric (diffractive) 
scattering a large number of trajectories representing 
pseudo-paths results. The Dyson integral equation can 
be converted to an algebraic matrix equation which 
can be solved by power series expansion or inversion. 
Using the examples of a circular and a stadium billiard 
we have shown numerically that the path-length power 
spectra calculated by the PSCA overcome shortcomings 
of standard semiclassical approximations by including an 
exponentially proliferating number of pseudo-paths and 
converges towards quantum transport. Moreover, we 
have presented a comparison between microwave billiard 
experiments [39], the SCA and the PSCA calculation 
and find evidence for contributions by pseudo-paths 
of PSCA to be present in the experimental data. We 
expect that our pseudo-path semiclassical approximation 
will be able to address unresolved issues of semiclas- 
sical ballistic transport such as the problem of weak 
localization [14, 34], the breakdown of symmetry of the 
autocorrelation function in reflection and transmission, 
and the semiclassical description of quantum shot noise. 



FIG. 9: (a) Test of unitarity U(k) on different levels of semi- 
classical approximations [SCA and PSCA (K max = 10)] to 
the rectangular billiard (D = L — \[A + n and d =0.25). The 
staircase function shows the QM result where U(k) = N(k) 
and N(k) is the number of open modes. Fig. (b) shows the 
discrepancy between the numerical integration of / and the 
FDA corresponding to Eq. (3.12) (solid line) and Fig. (c) the 
discrepancy in the phase angle. 



APPENDIX: Fraunhofer approximation and its 
limitations 

The integrals along the transverse coordinate (y) across 
the opening of the point contact are evaluated in Fraun- 
hofer diffraction approximation. The dependence of the 
action in zero magnetic field, kL a (y), on the transverse 
coordinate y is taken into account to first order Taylor 
series expansion, 



t /x ■ /„ s cos(6L.) 2 

L ai (y) = L ai + sm(6 at )y + y * y 2 ■ 



(Al) 



where L a = L a (y = 0) is the length of the path that 
starts from the center of the lead mouth. Keeping all 
three terms in (Al) leads to Fresnel diffraction integrals. 
Dropping the third term results in a Fraunhofer diffrac- 
tion approximation with the fundamental integral 
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d/2 



tFDA: 



,n) = 



2d 



a i{k n -\-k sin 0)y 



(A2) 



-d/2 



which, unlike Fresnel integrals, can be readily solved 
analytically in terms of elementary functions 



tFDA 



(9,n) 




ship; 



ksin9)d/2}\ 
- k sin 6 J 



(A3) 
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In terms of (A3) the reflection amplitude [Eq. (3.12)] is 
given by 

r(9, 6',k)= I FDA (6, 9', )k^2d(cos 9 cos 9') , (A4) 

with 9' = sin _1 (fc n /fc), the transmission amplitude 
[Eq. (4.6)] by 

t rn (9, k) = Vk^9 (I FDA {6, to) + I FDA (9, -to)) . (AS) 

The validity of the FDA hinges on the condition that 
the third term in Eq. (Al) is negligible. This term is of 
the order (kd)(d/L a ) < kd(jj) for short paths. Since in 
quantum billiards the asymptotic far-field regime is never 
reached, the FDA, somewhat counterintuitively, will fail 
for large kd. Indeed, in Fig. 9a we show a test of uni- 
tarity U(k) — R(k) + T(k) on different levels of semiclas- 
sical approximations to the rectangular billiard (inset of 
Fig. 9a) where the breakdown of the FDA and conse- 
quently of the PSCA for high modes to (m > 10) (or 
kd > 30) can clearly be seen. To highlight the failure of 
the approximation used in Eq. (Al) or (3.12) we plotted 



in Figs. 9b and 9c (solid lines) the discrepancy between 
the exact solution of the integral / [LHS of Eq. (3.12)] 
and the corresponding Fraunhofcr approximation. For 
L ai = L a2 = L and 9 ai = 9 a . 2 = as a function of 
k. Fig. 9b shows the difference in the absolute values 
|j| _ |jFZM| an( j Yig, 9c the absolute discrepancy in the 
arguments | arg(/) — axg(I FDA )\. Since the phase discrep- 
ancy reaches a fraction of unity for to > 16 (see Fig. 9c) 
random phases in the scattering amplitudes destroy path 
shadowing by pseudo-paths and therefore cause the vio- 
lation of unitarity (Fig. 9a). The failure of the FDA has 
more dramatic consequences for the PSCA than for the 
SCA. This is due to the fact that in PSCA the pseudo- 
path coupling leads to an exponentially growing number 
of FDAs with an increasing number of diffractive scat- 
terings, in contrast to the linear scaling of the number of 
FDAs involved within SCA. The point to be stressed is 
that this failure is not due to the PSCA but due to the 
additional FDA. Applying more accurate diffraction in- 
tegrals in this regime is expected to remedy the problem. 
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